To appear in Ap.J. 

The Micro- Arcsecond Scintillation-Induced Variability (MASIV) Survey III. 
Optical Identifications and New Redshifts 

Tapio Pursimo 

Nordic Optical Telescope, Apartado J^IJ^, 38700 Santa Cruz de La Palma, Spain 

tpursimoSnot . iac . es 
Roopesh Ojha^ 

NVI Inc./ U. S. Naval Observatory, 3450 Massachusetts Ave NW, Washington DC, USA 

David L. Jauncey 

CSIRO Astronomy and Space Science and Mount Stromlo Observatory, Canberra ACT 0200 

Australia 

Barney J. Rickett 

Department of Electrical and Computer Engineering, University of California, San Diego, La 

Jolla, CA 92093, USA 

Michael S. Dutka 

The Catholic University of America, 620 Michigan Ave., N.E., Washington DC 20064, U.S.A. 

Jun Yi Koay^ 

ICRAR, Curtin University, Bentley, WA 6845, Australia 
James E. J. Lovell 

School of Mathematics and Physics, University of Tasmania, TAS 7001, Australia 

Hayley E. Bignall 
ICRAR, Curtin University, Bentley, WA 6845, Australia 
Lucyna Kedziora-Chudczer 

^Now at NASA/Goddard Space Flight Center/ORAU, Code 661, Astroparticle Physics Lab., Greenbelt, MD 
20771, U.S.A. 

^Now at Faculty of Engineering and Science, Universiti Tunku Abdul Rahman, Setapak, KL 53300, Malaysia 



- 2 - 



School of Physics and Astrophysics, UNSW, Sydney NSW 2052 Australia 

and 

Jean-Pierre Macquart 
ICRAR, Curtin University, Bentley, WA 6845, Australia 

ABSTRACT 

Intraday variability (IDV) of the radio emission from active galactic nuclei is now 
known to be predominantly due to interstellar scintillation (ISS). The MASIV (The 
Micro- Arcsecond Scintillation- Induced Variability) survey of 443 flat spectrum sources 
revealed that the IDV is related to the radio flux density and redshift. A study of the 
physical properties of these sources has been severely handicapped by the absence of 
reliable redshift measurements for many of these objects. This paper presents 79 new 
redshifts and a critical evaluation of 233 redshifts obtained from the literature. We 
classify spectroscopic identifications based on emission line properties, finding that 78% 
of the sources have broad emission lines and are mainly FSRQs. About 16% are weak 
lined objects, chiefly BL Lacs, and the remaining 6% are narrow line objects. The gross 
properties (redshift, spectroscopic class) of the MASIV sample are similar to those of 
other blazar surveys. However, the extreme compactness implied by ISS favors FSRQs 
and BL Lacs in the MASIV sample as these are the most compact object classes. We 
confirm that the level of IDV depends on the 5 GHz flux density for all optical spectral 
types. We find that BL Lac objects tend to be more variable than broad line quasars. 
The level of ISS decreases substantially above a redshift of about two. The decrease is 
found to be generally consistent with ISS expected for beamed emission from a jet that 
is limited to a fixed maximum brightness temperature in the source rest frame. 

Subject headings: galaxies: distances and redshifts — galaxies: intergalactic medium — 
galaxies: ISM — galaxies: BL Lacertae objects : general — galaxies: quasars: general 



1. INTRODUCTION 



The discovery of centimeter wavelen gth, short terin intra-day variability in some compact, 
flat-spectrum extragalactic radio sources (|Heeschenlll984l : iHeeschen et al.lll987l ) was of immediate 
and profound as trophysical consequence . This phenomenon, which was soon dubbed Intraday 
Variab ility (IDV; Wagner &: Witzel 19951 ). implied brightness temperatures up to 10^^ K in extreme 
cases (jKedziora-Chudczer et al.l 119971 ) if the variations were intrinsic. Such enormous brightness 
temperatures present a serious challenge to the current paradigm for the physics of extragalactic 
radio sources. 
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However, there is now overwhelming evidence that IDV results primarily from scintillation 
in the turbulent ionized interstellar medium of our Gal axy, commonly refer r ed to as interstellar 
scintillation This explanation was proposed by iHeeschen &: RickettI (jl987l ) and has now 

been confirmed from several further lines of observational evidence. Time delays have been de- 
tected bet ween the arrival times of the intensity fluctuations from ISS sources at two wid ely-spaced 
telescopes (jjauncey et alJ l200d : IPennett-Thorpe de Bruynl I2OO2I : iBignall et alJ I2OO6I') . Further , 



annual cycles have been detected in the timescale of ISS sources (e.g. IJauncev &: Macquart 



2001 



Rickett et allboOll : iBignah et allboosl : IPennett-Thoroe k de Bruvnibood ). a periodic modulation 



that most plausibly results from the change in relative velocity of the Earth and the scattering 
medium through the course of a year. Subsequently, a strong correlation has been fo und between 



the level of IDV and Galactic emission measure estimated from Ha measurements (Lovell et al 



20081 ). This shows ISS to be the predominant c ause of radio IDV. Though some questions re- 



main ( Jauncev et al 



2001 



Krichbaum et al 



explanation of these observations (jJauncev et al 



20021). int erstellar scintillation is the only reasonable 



2001 



Though an extrinsic (as opposed to intrinsic) origin for IDV points to less extreme physical con- 
ditions for the sources that exhibit this phenomenon, scintillators are still among the most extreme 
and active radio AGN know n. For a source to scintillate, its angular size must be comparable to 
that of the first Fresnel zone (jNaravanlll992l ) which implies microarcsecond angular sizes for screen 
distances of tens to hundreds of parsecs. Such a high resolution cannot be achieved by any other 
existing technique, including space VLBI. Thus, inte rstellar scintillation has the potential to probe 
within a few light months of the central black hole (IBignall et al.l l2003l). Further, inferred bright- 
ness temperatures of a few scintillators are well in excess of 10^^ K (|Macquart et ahlboOOl ) which 



implies Doppler factors of several hundred or more (jReadhead 



19941 ) and is significantly higher than 



seen in VLBI surveys (e.g. iLister et al.ll2009l : lOjha et al.ll2010l ). While Doppler factors very similar 
to those found in VLBI are implied fr om the ISS analysis of the Green Bank radio monitoring 
program of 146 compact radio sources (jRickett et al.ll2006l ). an investigation of the properties of 
AGN that exhibit ISS is of considerable astrophysical interest, both because their extreme physical 
properties make them ideal probes of existing models of AGN physics and they probe phenomena 
at unprecedented resolution. 

The Micro- Arcsecon d Scintillation-Induced Variability (MASIV) survey was undertaken by 
Lovell et al.l (l2003l . l2008l ) in order to construct a large sample of scintillating extragalactic sources. 
A key discovery from the MASIV survey is that both the number o f sources showing ISS and their 
level of ISS variations appear to decrease with increasing redshift (jLovell et al.ll2008l ). This result 
was obtained for the 275 (of 443) sources with measured redshifts. While 206 of these redshifts 
were obtained from the published literature, 69 of them were based on preliminary analysis of our 
unpublished optical observations. A primary goal of the present paper is to publish the full analysis 
of those optical observations and an additional 22 sources subsequently observed. 



^In this paper we call the phenomenon of rapid variable flux density IDV and its cause, scintillation, ISS. 
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The apparent decrease in ISS with redshift is a new cosmologicahy important result. Its 
interpretation must include the cosmological prediction that the angular size of an object should 
initially decrease with redshift but start to increase above redshifts of about one. Observational 
evidence for this phenomenon has long been sought, but the best evidence to date is still only 



marginal (jGurvits et al.lll999l l. The MASIV sources were selected for their flat radio spectra and 
many are quasars. Thus any detailed interpretation must also address possible evolution in the 
beamed AGN jets presumably responsible for their very compact diameters. Alternatively, there 
could be scattering effects in the inter-galactic medium that could cause radiation from the more 
distant AGN to be scatter broadened so that their angular sizes are too great for them to scintillate 
in the interstellar medium of our Galaxy. 

A follow-up to the MASIV Survey has been conducted in which the ISS of a subsample of 140 
MASIV sources (70 with z < 2 and another 70 with z > 2) were observed sim ultaneously at 5 GHz 



and 8.4 GHz over a duration of 11 days using the VLA (jKoav et al.ll201ll ). These observations 
provided a means of determining the origin of the redshift dependence of ISS, by examining how 
the effect scales with frequency. This exploited the fact that intrinsic source size effects and scatter 
broadening both have different frequency dependences. The analyses and results were presented in 



Koav et al. 



torn . 



Our investigation of scintillators has been greatly handicapped by the lack of source redshifts. 
The MASIV survey found more scintillators among the weaker radio source sample, which also 
showed more extreme scintillation behaviour than their higher flux density counterparts. Since 
redshifts have been measured predominantly for the stronger radio sources, we particularly need to 
measure redshifts for the weaker sources. Redshifts are essential to determine physical properties 
including linear sizes, accurate brightness temperatures and luminosity distributions. Optical lumi- 
nosities are also an important probe of the physics of Compton scattering. For high Lorentz factor 
jets, energy losses should be dominated by Compton scattering of the diffuse radi ation field of the 



cosmic microwave background, starlight and reprocessed emission from the nucleus (jBegelman et al 



19941 ). This energy will escape mainly at X-ray wavelengths but some is expected in the optical 
band. Thus redshifts are needed to calculate the luminosities in order to understand the total 
energy budgets of these extreme sources. 

In this paper we present new redshifts obtained with the 2.56 m Nordic Optical Telescope 
(NOT) on La PalmsEI and the 5 m Hale Telescope at Mount Palomar. We also present redshifts 
obtained from the literature whose reliability could be ascertained and optical magnitudes collated 
from the surveys and literature. Further, we look for relationships between the redshift, spectral 
classification, luminosity and the observed level of ISS variation and compare the results with a 
simple model for ISS. 



^Based on observations made with the Nordic Optical Telescope, operated on the island of La Palma jointly by 
Denmark, Finland, Iceland, Norway, and Sweden, in the Spanish Observatorio del Roque de los Muchachos of the 
Instituto de Astrofisica de Canarias. 
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Throughout this paper we adopt a cosmology with = 0.27, = 0.73 and Hq = 70 km s 
Mpc~^. We use the convention S,^ oc for flux density at frequency v. 



2. THE MASIV SAMPLE 



The MASIV sample has its roots in two radio catalogues: CLASS (C osmic Lens All Sky 
Survey; Myers et al . 2003) and JVAS (Jod rell Bank VLA Astrometric Survev: IPatnaik et al.lll992l : 
Browne et al.lll998l : IWilkinson et al.lll998l ). Both of these were targeted surveys of flat-spectrum 
radio sources using the VLA at 8.4 GHz in the A' configuration. JVAS is eff'ectively a bright sub- 
sample of the complete J VAS/CLASS sampl e, which consists of more than 11,000 sources with flux 
density down to 30 mJy (jMvers et al.ll2003l ). We note that the CLASS flux density limit is well 
below the 100 mJy flux density limit chosen for the MASIV sample (see below). 

For MASIV observations, we wanted to target compact flat-spectrum sources that were un- 
resolved at 8.4 GHz with the VLA and located north of the equator. As a first step, all CLASS 
sources with modeled source sizes < 50 mas and all JVAS sources with > 95% of their fiux density 
in an unresolved component were selected. These selection criteria were chosen due to the way the 
data were presented in each catalogue: for JVAS, the peak and total VLA flux densities are listed; 
for CLASS the modeled component size is presented. Both criteria select those sources which are 
essentially unresolved point sources for short (0.5 - 1.5 minute) snapshot observations with the VLA 
in A' conflguration (this is the highest resolution configuration of the VLA). Sources where the 
catalogues indicated the presence of any confusing sources - indeed, any emission outside of the un- 
resolved target source - in the VLA field, were dropped from the sample, as selecting each MASIV 
target field to consist of a single, isolated point source greatly simplifies the analysis of variability. 
To obtain a flat-spectrum sample with spectral index a > —0.3 ( S,y oc z^"), these so urces were 
cross-correlated with the NVSS catalogue (NRAO VLA Sky Survey: ICondon et al.lll998l ). although 
we note that the observed flux densities in these catalogues were non-simultaneous. 

The sample was next divided into strong and weak subsamples of about 300 sources each with 
the strong sample consisting of sources with Ssaghz > 0.6 Jy and the weak sample consisting of 
sources with 0.103 < Ssaghz < 0.13 Jy. Finally, the sources were selected to have uniform RA - 
6 distribution in o rder to have the best possible coverage for the VLA observations (for details see 
LoveU et al.1l2003l ). 



MASIV's 5 GHz observations were carried out on this 'core' sample of 578 sources. Of these, 
102 sources were removed in the analysis stage due to the presence of structure or confusion while 
one was removed due to an error in the initial sample selection process. This left a sample of 
475 point sources common to all four epochs of MASIV observations. However, as described in 
Lovell et al.l (120081 ). 32 of these sources were used as secondary calibrators in two or more epochs, 
which effectively removes their low- level (instrumental or real) variations in those epochs, relative 
to the other sources. These 32 sources were excluded from the final MASIV sample, leaving 443 
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sources for full analysis. 



As the average flux densities of many sources from the originally defined radio "strong" and 
"weak" sai nples vary with re spect to the original catalogued flux densities, we used the first epoch 
VLA data (iLovell et al.ll2003l ) to set the dividing flux density to 0.3 Jy at 5 GHz. In the rest of the 
paper we refer to the radio strong sample of 229 sources (including the 32 secondary calibration 
sources for the optical analysis only, Section 4) and the radio weak sample of 246 sources. 



Lovell et al.l (|2008l ) report which of these sources exhibited Intra Day Variability (IDV) and 



in how many epochs of observation were classified as variable. A structure function (SF) analysis 
of the flux density variations was also done, which provides a quantitative measure of the strength 
of the variation. The flux density was normalized by its mean and the average SF was computed 
from the 4 epochs and tabulated for each source as I?(2days), the SF at a time lag of two days. 
This quantifies the variations and was shown to be correlated with the interstellar electron column 
density, estimated from Ha emission on that line of sight. The correlation confirms that the 
variations are indeed predominantly due to ISS. There are thus two ways of quantifying the ISS 
of a source, either as the number of epochs in which it was classified as variable or the 4-epoch 
mean D(2days), where w e note that no sigii ificance should be attached to values smaller than about 
4 X 10"'' (for details, see lLoveh etliLlbnOsI ). 



3. REVIEW OF THE LITERATURE 

The main goal of this work is to present the spectroscopic identifications and redshifts of the 
flat-spectrum, extragalactic radio sources that make up the MASIV sample. 



3.1. Redshifts from the Literature 

We searched for redshifts and spectroscopic identifications of MASIV objects in the litera- 
ture starting from data repositories like the NASA Extragalactic Database (neeI). Spectroscopic 
identifications and redshifts were found for the majority (171 of 229) of the radio strong sample 
sources. However, this is not the case with the weak sample (62 of 246). Most li terature red- 
shifts and spectroscopic identifications are from the 12th Veron- Getty Sz Veron (j2006l ) compilation 
of known AGN, hereafter yGV12. Identificat ions for 27 objects are from GGRaBS (Gandidate 
Gamma-Ray Blazar Survey: iHealey et al.ll2008l ) and the rest from sources found using NED. 



For those objects which had reported redshifts, we tried to locate their spectra and/or infor- 
mation about their emission lines (line identifications, equivalent widths) from which the redshift 
estimates were obtained. We classified the sources into categories A, B or G, based on the reliability 



• ^http: //nedwww. ipac . caltech. edu/| 
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of their redshift estimates. Sources with redshifts obtained from multiple reported emission lines or 
with spectra clearly showing multiple emission lines were placed in category A. If only a single line 
was reported with wide wavelength coverage, or the spectrum provided was of marginal quality, 
the source was placed in category B (see comments in Table [2l). Redshift estimates of sources for 
which the spectra were not available, but where the process by which the redshifts were obtained 
were clearly described, were also placed in category B, as were three sources where imaging studies 
have detected the host galaxy but show featureless spectra. For these three objects, we adopted the 
imaging redshift whic h is estimated assuming the host galaxy is a "standard candle" (for details see 



Sbarufatti et al.ll2005l ). Finally, objects in which the redshifts were indicated to be photometric, or 



where no further information could be found on how they were obtained, were placed in category 
C. Objects with (nearly) featureless spectra, with conflicting reported redshifts, or with redshifts 
estimated based on an intervening absorption system, were all placed in category C. Comments on 
archival redshifts that either were rejected or judged less reliable are in Appendix A. 



3.2. Optical identification 



In order to define the sample for optical spectroscopy and select appropriate observing re- 
sources, accurate optical identification was a prerequisite. 

For optical identification we used the CDS client databas^ to access three large area surveys, 
the Sloan Digital Sky Survey (SPSS), the Guide S tar Catalog 2.3 (GSC) and the USNO-Bl Catalog. 
The SDSS DR5 (jAdelman-McCarthv et al.ll2007l ) is the primary catalogue for optical ide ntification. 



It is complete down to r=22.5 (Ag// ~6230A) with astrometric accuracy of about (/.'l ( Pier et al 
20031 ) and covers essentially the Galactic gap (RA ~12'^ 6 ~30°) and some smalle r patches altogether 



about 8000 square degrees. As a secondary catalogue we used the GSC 2.3 (ILasker et al.l l2008l ) 
which is compiled from scans of the POSS-II plates and has F magnitudes (Ag// ~6500A) down to 
~20.5. This has poorer photometric (0.13-0.22 mag) and astrometric accuracy (0'.'3-0'.'4) than the 
SDSS. For 16 sour ces which were not found from the GSC2.3 we used USNO- Bl R2-magnitude s 
(jMonet et al.ll2003l ). which has poorer photometric accuracy than the GSC2.3 (jSesar et al.ll2006l ). 
Of these sources ten are faint {R2 >19.0), just visible from the POSS-I 0-plates, and four are 
bright extended targets {R2 <11.3). The ten faint objects might be missing from the GSC2.3 due 
to variability, while the four bright, extended objects may be missing due to the design of the 
GSC2.3. Finally, for eight faint sources we found optical magnitudes from the literature. 

The initial search radius around the radio position was 5". From SDSS DR5 we found 165 
matches to a MASIV radio source within one arcsecond, with median r-magnitude of 18.7ibl.8 
(standard deviation of the distribution) and the faintest optical identifications having r=22.9. The 
mean separation increases from O'.'O? for r < 20.5 mag sources to O'.'IS for the 37 sources with 20.5 



' jhttp: //cds .u-strasbg. f r/ 
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< r < 22.9. The GSC2.3 has 356 matches within I'.'S with median magnitude of 18.4ibl.5. The 
mean separation of F < 19.5 sources is 0'.'25 and for the 38 sources with 19.5 < F < 20.7 it is 
(y.'39. Fewer than 10% of the GSC2.3 matched sources have separations greater than 0.5 arcsecond 
and only ~2% have separations greater than one second of arc. In Figure [1] we show a plot of 
the measured separation versus apparent magnitude. In both the SDSS DR5 and GSC2.3 matched 
sources, the accuracy of the astrometry is as expected and decreases close to the limiting magnitude 
of the survey. 

Multiple matches were found in 14 cases (~4% of the matches). Inspection of the data and 
visual inspection of the images (available from the CDS) revealed that close companions (8 sources), 
defined as those with a separation less than one arcsecond, were from separate epochs. This suggests 
that there is only a single source, but the astrometric solution is slightly different from epoch to 
epoch. For the remaining six sources the nearer object was selected as the optical counterpart and 
a second target was found with the angular distance 2.5-4.9 arcseconds from the radio source. Of 
the 356 GSC2.3 detected objects, 133 also have SDSS DR5 magnitudes, hence we are left with 223 
sources with only F-magnitudes. Down to a Galactic latitude of 5 degrees, all SDSS and GSC2.3 
matches had only a single source within the search radius. Table [J summarizes the SDSS and GSC 
identifications. 

In addition we defined a subsample, (08'^ <RA< 16^ and 5 < 64°), which covers roughly 
the SDSS DR5 sky area (hereafter 'SDSS sample'). Of the 79 radio strong sources, all but one 
secondary calibrator source have optical identifications. Optical identifications are available for 82 
out of 95 sources in the radio weak sample. 



3.2.1. Combining GSC2.3 and SDSS DR5 Data 



We transformed the GSC2.3 and SDSS DR5 data to the common magnitude system using 



RsDSS 



0.272(r-i)- 0.159 



(1) 



bv IChonis fc Gaskelll (120081 ) for the SDSS data. This is derived using SDSS DR5 and iLandolt 
(|l992l ;) data for r > 14 magnitude stars. For the r < 14 stars they found a systematic magnitude 



difference, however we have only two such bright targets. 

We obtained the G SC2.3 F- and N -magnitudes for the standard stars, which were observed 
more than four times by iLandoltl ()2009l ) where we found (almost) a one to one correspondence 



Rgsc 



(2) 



The F — N color term and zero point corrections are small and are omitted (see Appendix B for 
details). The USNO-Bl R2-magnitudes were converted in a similar way as the F magnitudes. For 
249 sources with 13 < < 19 the average magnitude difference F — R2 is ~-0.07. Taking into 
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account the small number of -R2-magnitudes and the magnitude offset, the possible bias to the 
global properties is negligible. 

We checked the common magnitude system by selecting all targets with SDSS photometry and 
GSC2.3 F magnitude between 16.5 and 19.5, altogether 97 sources. The median r - F magnitude 
difference is 0.3ib0.4 and after the Rsdss ~ Rgsc transformations the median difference is O.lzbO.4 
(Figure [2]) suggesting that the systematic difference between the transformed Rsdss and Rose 
magnitudes is small in comparison to the uncontrollables, such as variability of the sources and 
photographic plate sensitivity changes. 

The magnitudes are integrated magnitudes and there has been no attempt to separate the host 
galaxy contribution from the nucleus or correct the emission line fluxes. For optically extended 
sources (all low redshift, z < 0.2 sources as well as Sy2/NLRG-type objects at higher redshift) the 
absolute magnitude includes both the core and the host galaxy. For example, J0057+3021 (NGC 
315) has USNO R raagnitu de ~11.1 however the core R-magnitude is ~19.8 (assuming V~R=0.3; 
Verdoes Kleiin et al.ll2002l ). Also, the GSC2.3 photometry pipeline was tuned for point so urces, 
thus there is a systematic overestimate for bright non-stellar magnitudes (jLasker et al.ll2008l ). 



All the magnitudes w e re de -reddened using the E{B — V) values and the relative extinction 
values from ISchlegel et aP ( 1998 ). In Table [2] and for the optical luminosities, we refer to Rsdss 
and Rgsc as i?-magnitude. 



3.2.2. Radio and Optical Properties 

The optical identification rate, combining the results from the SDSS DR5 and GSC2.3 cata- 
logues, is ~ 86% for the radio strong and ~ 77% for the radio weak sample. However, since the 
two catalogues have different sky coverage and hence different Galactic latitude coverage, this could 
bias the identification rate. Therefore, using only the (all sky) GSC2.3 catalogue, the strong sample 
has ~ 80% (187) of F-magnitudes brighter than 20 compared to ~ 60% (146) for the weak sample. 

Including empty fieldjf], the median optical flux density of the strong radio sample is ~0.15 
mJy (18.2 mag) decreasing by a factor of 3.2 to ~0.047 mJy (19.5 mag) for the weak sample. In 
comparison, the median radio flux density decreases by a factor of 7.4, from ~0.89 Jy to ~0.12 
Jy in the two samples. If we were to use this radio flux ratio to calculate the weak radio sample 
median optical magnitude from the strong radio sample mean optical magnitude, it would yield a 
i?-magnitude of 20.1 for the weak sample. The /c-correction difference between the weak and stron g 



samples is unknown, but most likely it is ~0.2 magnitudes or less on average (e.g. IWisotzkil l200d ) . 
assuming the average redshifts of 1.2 and 1.8 for the strong and weak samples respectively (see 
below) . 



^In the optical there is a well defined upper detection limit, i.e. the median is a non-biased estimator. 



- 10 - 



Figure [3] shows the optical versus radio flux density with diagonal lines indicating constant 
radio-optical spectral index (defined by Sy oc 1/°'^°), which has only been corrected for Galactic 
extinction and lacks any /c-correction because 40% of the radio weak sample objects have un- 
known redshift. Th e top line corresponds to 0^0=— 0.2, the generally accepted radio quietness limit 



(jStocke et al.lll985l ). Although there appears to be no detailed correlation between optical and radio 
flux density, the weak sample is on average one optical magnitude fainter than the strong sample. 
The radio weak sample also appears to have an even larger spread than does the radio strong 
sample. The weak sample has 17 sources with a^o > —0.45 and the strong one only five. Close 
to aro=-1.0 both samples have many optically empty fields thus we cannot draw any conclusions 
about the underlying distribution. It is also important to note that the non-simultaneous nature 
of these radio and optical observations will further increase the scatter. The difference between the 
radio to optical flux ratios of the radio strong and weak samples is statistically significant. Running 



the ASURV Kaplan-Meier estimator (|Avni et al.lll980l ). the probability that the two samples are 



drawn from the same parent population is small (P <0.01%). This holds for the 'SDSS sample' 
and also for the whole sample using a shallower GSC2.3 limiting magnitude. 

Our results indicate that the radio faint objects are optically "brighter" than expected using 
the radio strong objects as reference. The explanation for the difference is not clear. One possibility 
is that the radio weak sources have larger host galaxy contributions in optical than radio strong 
objects. Another possibility is that the average spectral energy distributions of the strong and 
weak sources are different, perhaps due to differing contributions of the accretion disk ("big blue 
bump") to the optical flux density or different amounts of dust. 



4. SPECTROSCOPIC FOLLOW-UP 

The new redshifts reported in this paper were obtained using the 2.56 m Nordic Optical Tele- 
scope (NOT@) which is located at Roque de los Muchachos, La Palma, Canary Islands, Spain and 
the 5 m Hale Telescope at Mount Palomar, California, USA. Most of the NOT data were obtained 
during two observing runs (July 3rd to 7th, 2005 and July 22nd to 28th, 2006) supplemented by 
a few additional nights between January 2004 and July 2007. The Mount Palomar data were ob- 
tained during an observing run from August 9th to 10th, 2007. Most fields were pre-imaged at 
NOT with deeper and better image quality than the SDSS frames to verify the optical identification 
of the radio source. Altogether 90 objects were observed at NOT and 28 objects were observed at 
Palomar with some targets observed at both telescopes. 

At NOT, low resolution spectra were obtained using the Andalucia Faint Object Spectrograph 
and Camera (ALFOScjll) with Grism 4 (300 rulings/mm, 3200-9100 A, dispersion of 3A/pixel) with 



|http: //www.not . iac . es 

The data presented here have been taken using ALFOSC, which is owned by the Institute de Astrofisica de 
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second order blocking filter giving wavelength coverage about 3800-9100 A. However, in practice the 
NOT red limit of the spectrum is reduced to about 8000 A due to strong fringing. We used Grism 
6 (600 rulings/mm, 3200-5550 A) or Grism 3 (400 rulings/mm, 3200-6700 A) for objects with a 
single detected line at about 5000 A in order to confirm the redshifts from either Mgll, GUI] or 
the Ly-alpha line. The NOT run in 2005 was photometric, based on Garlsberg Meridian Telescope 
(GMT) datsl^, with seeing of 0'.'6-l'.'4 and the 2006 run had some dusty and/or non-photometric 
nights with O'.'7-l'.'l seeing. At Hale, the Double Spectrograph (DBSP) was used with a setup of 
600 lines on the blue arm (3400-5400A, ~lA per pixel) and 158 hues for the red (5400-9500A, 
~4.9A per pixel) with either I'.'O or 1'.'5 slit giving continuous wavelength coverage from about 
3400 to 9500 A. 

Most of the targets were observed with the slit in parallactic angle in order to reduce the light 
losses in blue. We note that, each spectroscopic target has a characteristic AGN spectrum and there 
is no evidence of any spurious counterparts. In two cases the nearby companion was put on the slit 
and in both cases the object closer to the radio source had a characteristic AGN spectrum while 
the companion had a typical stellar spectrum. For the NOT data, an internal flat field (Halogen) 
image was obtained before each science exposure, in order to improve the fringe correction at the red 
end. In addition, several HeNe images were taken during the night for the wavelength calibration. 
The DBSP calibration frames, dome flats, HeNeAr (red arm), FeAr (blue arm) images, were taken 
morning and evening. The data were reduced using standard IRAF§ procedures including bias 
subtraction, flat field correction using internal flat field lamp and wavelength calibration. Before 
extracting the spectrum, cosmic rays were removed from individual spectrograms. The majority of 
the spectroscopy targets have R-magnitude between 17 and 21, with a median of 18.8. 

The redshifts were determined from the narrow emission lines (e.g. [O IIlJAA 4959, 5007, [O II]A 
3727) whenever possible, or from strong broad emission lines otherwise, typically H /3, Mgll , G III], 
G IV or Lya. The new redshifts are based on two or more line identifications. However, in four cases 
with wide spectral coverage and high S/N, we assumed the line to be Mgll , as if the line were to be 
G III], we would have expected to detect G IV] or Mgll . In Table[2]we list, respectively, source name, 
RA, declination, first epoch 5 GHz fiux density, epochs of variability (at 5 GHz), Z)(2days), Galactic 
reddening, de-reddened R-magnitude, R-magnitude reference, redshift, spectroscopic identification 
with the original references and a comment. Selected objects are discussed in Appendix A and G. 
Objects with only a single detected emission line are listed in Table El 



Andalucia (lAA) and operated at the Nordic Optical Telescope under agreement between lAA and the NBIfAFG of 
the Astronomical Observatory of Copenhagen. 

: //www. ast . cam. ac .uk/$\sim$dwe/SRF/camc$_$extinct ion. html | 

^IRAF (Image Reduction and Analysis Facility) is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in Astronomy (AURA) under cooperative agreement 
with the National Science Foundation (NSF). 
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4.1. Classification Scheme 



On the basis of the optical spectrum the sources have been divided into three main groups: 
Type objects with weak emission lines (mainly BL Lacs), Type 1 objects with strong and broad 
emission lines (mainly FSRQ) and Type 2 objects with strong and narrow emission lines (mainly 
Seyfert 2's). In the context of AGN unification schemes, Type and Type 1 objects have their 
Doppler boosted relativistic jet aligned close to the line of sight with Type 2 objects seen at a 
larger angle in such a way that the broad line region is obscured by a dusty torus (see review by 
Urrv &: Padovanilll995l). The spe ct roscopic classifica ti on sche me is adapted from pre v ious w orks e.g. 
Laurent-Muehleisen et al.l (119981 ) , ICaccianiga et al.l (l2002bl ) and ICaccianiga et al.l (120081 ) . Many 
observed properties of the AGN, such as emission line equivalent widths (EW) and luminosities 
have continuous distributions. Hence the division between, for instance. Type and Type 1 based 
on equivalent widths EW or between Seyfert 1 and FSRQ based on optical luminosity depend on 
the dividing line chosen. The VCV12 classifies all objects fainter than Mb = —23 (using VCV12 
cosmology and magnitudes) as (low luminosity) AGN, which is converted to our cosmology and 
R-band resulting in a QSO vs radiogalaxy (RG) limit Mr = —22.8. As noted by the VCV12, the 
Mb =-23 is an arbitrary limit without any physical reason. Due to different source classification 
schemes, and/or to variability, a QSO might have been classified as an AGN in one epoch but as 
an FSRQ at another e poch, and vice versa . For exam ple, a strong radio source, J 1 506+3730 is 
classified as an AGN bv lHealev elTall NLRG bv ISowards-Emmerd et aD ^2Qo4 ) and Sy2 by 

VCV12. Also, in some cases the classification depends on S/N, e.g. a weak AGN-core wit h dilution 
by th e host galaxy where broad emission lines might be undetected; f or J 0057+3021, 



Ho et al 



(119951 ) reported a broad Ha line, but a later study could not confirm it (jGongalves fc Serote Roos 
20041 ). 



We describe our classification scheme below. Note that all the line equivalent widths (EW) 
are rest frame values: 



Oa) BL Lacs: objects with no lines or line rest frame EW < SA and C < 25%, where contrast C 
is defined as 

C = {U- /_)//+, (3) 

where /_ is the mean flux in the 3750-3950 A range and /+ in the 4050-4250A range where 
flux is measured in frequency. In the case of moderate S/N we give only tentative spectroscopic 
identiflcation. 

Ob) BL Lac candidates: 25 < C < 40% and strongest line EW < 40 A. 
Oc) Passive elliptical galaxies (PEG): C > 40% and EW < 40A. 

la) Flat spectrum radio quasar (FSRQ): at least one permitted line excluding Ha, line FWHM > 
1000 km/s and Mr < -22.8. 

lb) Seyfer t type 1 galaxie s: Mr > -22.8, the subclass has been defined based on the Hf3/[0 III]A5007 
line ratio (Iwinklerlll992 ). 
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Ic) Narrowline Seyfertl/QSO (NLSyl/NLQSO): At low redshift: Hp/[0 III]A5007 < 3 and Hp < 
2000 km/s. For the z > 1.5 targets: rest frame line width 1000-2000 km/s and th e relative strength 



of the He II emission line compared to C IV following iHeckman et alj (|l995l ): ICaccianiga et al 
(|2008l ). Type 1 objects have weak He II emission line. 



2a) Narrow line Radio galaxy (NLRG): line FWHM < 1000 km/s and Mr < -22.8. 

2b) Seyfer t type 2 galaxie s: Mji > -22.8, the subclass has been defined based on the H(3/[0 III]A5007 

line ratio (jwinklerlllQOa ). 

2c) Galaxy: C > 40%, absorption lines e.g. Mg I and no emission lines 

Low redshift sources are the most problematic to classify as host galaxy dilution can hide 
broad lines and can reduce the line EW. Also the non-thermal continuum can be difficult to detect. 
Fortunately only four of our follow-up targets are at low redshift {z < 0.2), suggesting that the bias 
from the uncertain spectroscopic identification is small. 



4.2. Summary of the spectroscopic follow-up 

We report 79 new redshifts with spectroscopic identifications, 7 new BL Lacs with featureless 
spectra, and twelve objects with a single emission line of which six have a tentative redshift (Table [3] 
and Appendix C). We confirm the BLLac status of nine objects, and confirm the archival redshifts 
as well as spectroscopic identifications for nine other objects. We also re-classify two BL Lac objects 
using our own data (Appendix A). Spectroscopic identifications from the literature were reviewed 
using the same scheme as for our own data, whenever possible. In some cases this led us to adopt 
a different classification than NED (see appendix A) . Example spectra of our own observations are 
shown in Figured! 



4.3. RESULTS AND SAMPLE PROPERTIES 

Here we present 347 MASIV sources which have their spectra classified, either from the lit- 
erature or based on our own data (Table [2]). The radio strong sample is 91% (208/229) complete 
and the radio weak sample is 57% complete (139/246). Of the twenty-one radio strong sources 
without spectroscopic identification, seven have magnitudes from the surveys or literature, with 
R-magnitude between 17.3 and 23.99. 

In TableHlwe summarize the number of sources in each optical spectral classification, separating 
the radio strong and weak subsamples. It shows that Type 1 sources are the most common, 
comprising 78% of both the radio strong and weak samples. About 18% of the radio strong sources 
and 16% of the radio weak sources have Type spectra. A small minority of the sources have Type 
2 spectra, 4% and 6% of the radio strong and radio weak samples, respectively. It also shows that 
almost all Type 1 sources are optically luminous FSRQs. In Table [5] we compare the identifications 
based on the present sample and the 'SDSS sample'. The results suggest that the two samples have 
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similar distributions of the sources for both radio strong and weak subsamples. Also the redshift 
distributions are identical. This suggests that the observed sample is representative in comparison 
to the full 475 object sample. 

Though broadly similar, there is some difference in the spectroscopic identifications of the radio 
strong and weak samples. The radio weak sample appears to have slightly more low luminosity 
AGN (e.g. PEG, LINER, Hll-region, Galaxy) and fewer BL Lac objects than the radio strong 
sample. An increasing number of low luminosity objects are expected to be detected as the flux 
limit of the sample decreases. At low redshift {z < 0.2) most (7 out of 8) of the radio strong sample 
objects are BL Lacs or have broad emission lines, in contrast to the weak sample, where half (6 
out of 12) of the z < 0.2 objects have narrow emission lines or "galaxy" type spectra and a further 
five are classified as PEG. Also, the PEG and galaxy type objects are absent in the radio strong 
sample. 

All together, redshift information is available for 319 sources. Figure [5] shows the redshift 
distributions of subsamples divided by the optical spectral classification and the MASIV first epoch 
5 GHz flux density. For Type 1 sources the median redshift increases from 1.3 to 1.4 as the median 
(5 GHz) flux density decreases from 0.86 to 0.12 Jy. The Type 1 radio strong sample has a clear 
maximum in the distribution, but the radio weak sample has a fairly flat distribution between 0.7 
and 2.0. Of the Type 1 sources ~20% of the radio strong and ~29% for the radio weak sample 
have z > 2.0. The redshift trend due to the flux density limit is similar to that seen from the 
previous blazar surveys e.g. from S 5 Zm.pan= 1-18 (li mit IJy) to S4 Zmean= 1-3 (limit 0.5 Jy) to 



DXRBS Zmean=l-56 (fimit 0.05 Jy) (ILandt et al.ll200lh . Most of the Type or Type 2 sources have 



measured redshifts of less than one and almost 50% of the Type sources are without a redshift. 

We note that there might be a small number of mis-classified objects. As in every survey, the 
signal to noise ratio of the spectroscopic data is not always optimized for spectroscopic classifica- 
tion but for measuring the redshift. However, in most cases there is no ambiguity in our source 
classification. As our main results and conclusions are based on Type 1 sources (see below), the 
effect of some mis-classified objects is not significant. 



4.4. MASIV and other Blazar Surveys 



In Table [6] we compare spectroscopic identifications, limiting magnitu des and the radio fiu x 
density limit of MASIV with some other blazar surveys. The DXRBS sample ([Padovani et al.ll2007l ). 
i s selected by cross -correlating X-ray sources with radio sources {S^ghz ^ 50m Jy). The CGRaBS 
(iHealey et al.ll2008l ) has fiat spectrum radio sources selected from the CRATES (iHealey et al.ll2007l ) 
survey with S^^sGHz > 65 mJy, where the aim is to provide a catalogue of likely 7-ray loud AGN. 
This selection is based on a figure of merit (FoM) which measures the likelihood that an individual 
radio/X-ray source is associated with a kiiown 7-ray source. In term s of selection criteria the CLASS 
blazar survey (CBS) (jMarcha et al.ll200ll : ICaccianiga et al.ll2002bl ) is most similar to MASIV. The 
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CBS has flat radio spectrum sources with the weakest sources S^ghz ^ 30mJy and red magnitude 
< 17.5. It is interesting to ask how the source identification differs between the MASIV selection 
criteria in comparison to the previous blazar surveys. 

The DXRBS classification into the categories of "Radio Loud Quasar, BL Lac and Radio Galax- 
ies" is almost identical to our Type 1, Type Type 2. The possible difference is that they might clas - 



sify "PEG" as "Radio Galaxy" and their survey is nearly complete (^95%. iPadovani et al.ll2007l ). 
The CGRaBS classifies objects as FSRQ (Type 1), AGN (Type 1), BLL(Type 0), NLRG(Type 2) 
and GALAXY (Type 2) where our corresponding classifications are in parenthesis. CGRaBS is 
about 79% complete with respect to the entire survey a nd 85% for objects with kn own R-mag< 
23. The comparison between the CBS (~91% complete: ICaccianiga &: Marchal 12004 ) and MASIV 
is straightforward, as the MASIV classification is adopted from the CBS. 

The distribution of the spectroscopic classes is almost identical between MASIV and the 
DXRBS and the CGRaBS with slightly more broad line quasars in MASIV than other surveys. 
The main difference between the CBS and MASIV is that the CBS has more low redshift Type 
2 objects than MASIV. This difference remains when using similar radio ('CBS rb') and optical 
flux density ('MASIV ob') limits for these samples. This indicates that the selection of compact 
radio structure in deflning the MASIV sample, filters out nearby Type 2 objects and favors FSRQs. 
Probably as a result of its bright optical fiux density limit, the CBS has at least 40% of low redshift 
{z <0.15) sources. In contrast, MASIV has 4% for the total sample, 1% and 6% for the radio strong 
and weak samples, respectively. In all these cases it is assumed that the BLLacs without redshift 
have z >0.15. Another difference between MASIV and the CBS is that sources that are optically 
less luminous but powerful in the radio (Figure[6]top panel, where Mji > —23 and P > 10^^ W/Hz) 
are missing from the optically bright CBS sample. 

The source and the redshift distributions of the MASIV, DXRBS and CGRaBS samples are 
similar. There might be slightly fewer z > 3 sources in CGRaBS (~ 3%) than in MASIV (~5%). 
The difference could be, for example, due to incompleteness of the MASIV sample or pre-selection 
of only compact radio sources for the MASIV sa mple. MASIV has fewer low radio power FSRQ 



sources (P^ < lO^^-^ W/Hz) than DXRBS (22%: iLandt et al.ll20nih . but about the same amount 
(~2-3%) as the combined 1-Jy and S4 samples. Some of the sources classified as "galaxy" might 
have weak broad lines which are not detected due to dilution by the host galaxy continuum and/or 
insufficiently high signal-to-noise ratio of the data. 

Apart from decreasing the number of low radio power sources, using compact radio core selec- 
tion increases the number of objects in the MASIV sample that have featureless optical spectrum 
i.e. are classified as BL Lacs. Specifically, it is interesting to note that IDV4 (i.e. IDV detected in 
all four MASIV epochs) is one of the most effective preselections to find BL Lac objects. The use 
of persistent IDV as a selection criterion strongly favours BL Lac objects with 43% of the MASIV 
sources showing IDV in all four e pochs identified as BL L ac objects. In comparison, a radio- X-ray 
selection in the XB-REX sample (jCaccianiga et al.ll2002al ) and DXRBS (jPadovani et al.ll2007l ) have 
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a BL Lac fraction of 25% and 18% respectively. However, 7-ray loudness may b e an even stronge r 
selector of BL Lacs with about 50% of Fermi/LAT identifications being BL Lacs (jAbdo et al.ll2009l ). 



Lister et al.l (j200ll ) studied the IDV properties of the Pearson- Readhead compact extragalactic 



radio sources and found that IDV sources have smaller emission line widths and lower 5 GHz 
luminosites than non-IDV sources. Their results about line widths are consistent with our findings 
that 25 to 40% (for IDVl through IDV4 subsamples) of the MASIV IDV sources are BL Lac 
objects and that 70% of the BL Lacs are IDV sources. Whether there is a correlation between 
IDV and line width amongst the Type 1 sources will be studied in a future paper. The 5 GHz 
luminosity difference found by Lister et al. could be due to combining FSRQs and BL Lacs into one 
IDV/non-IDV sample, and the different redshift distributions of the IDV and non-IDV samples. 



5. Discussion 

5.1. Dependence of ISS on 5 GHz flux density and optical spectral type 

The correlatio n between IDV and the 5 GHz flux density w as discovered in the MASIV survey 



(jLovell et al.ll2003l ) and later confirmed by lLovell et al.l (120081 ). In the latter work, the variability 



of the flux density was studied using two complementary methods (see Section 2): D(2days) which 
combines all the data in a statistical estimator; and a visual analysis of each light-curve at each 
epoch, which classified each source at each epoch as variable or not variable and so allowed an 
estimation of the duty cycle in the variability. Sources were classified as "ISS" if they were found 
to be variable in 2 or more epochs and as "non-ISS" otherwise. 

To study the flux density/ISS correlation in detail, we divided the sample by the flux density 
and spectroscopic type. Figure [7] shows the -D(2days) histograms of the MASIV subsamples. There 
is a clear increase of the median and mean Z)(2days) as the flux density decreases and this is seen in 
all spectral types separately. One can see that many of the distributions in Figure [7] are significantly 
skewed, having a tail extending to high values. When binned on a linear scale the distributions 
are more strongly skewed. A measure of this skewness is that the mean values of D(2days) are 
considerably larger than the median, since the rare high values influence the mean but have little 
effect on the median. 

The Kruskal-Wallis (KW) test indicates a ~0.4% probability that the radio strong and radio 
weak Type samples have equal median Z)(2days). The probability that Type 1 radio strong and 
weak samples have similar median is negligible {pkw < 0.01%). Similarly, the Kolmogorov-Smirnov 
(KS) test indicates that the Type radio strong and weak distributions of D(2days) are drawn 
from the same parent population with a probability of 0.1%; for Type 1 sources this becomes pxs < 
0.01%. For Type 2 objects the KW and KS tests indicate the same median or parent population 
with less than 5% probability, however the full sample contains only 19 objects and their D(2days) 
values are typically at or below the threshold (0.0004) for significant variation. This indicates that 
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the -D(2days) has a S^ghz dependence and this holds for Type and Type 1 sources with high 
statistical confidence. A similar, but less significant, trend is also seen for Type 2 sources. We note 
that we have similar results when comparing the 5 GHz flux density distributions for sources drawn 
from the 'SDSS sample'. Thus we confirm the lYyV-S^cHz flux density dependence found earlier 
and show that this correlation is present at least for Type and Type 1 sources. 

In previous studies of IDV, the optical spectroscopic type of the objects has not been consid- 
ered. Figure [7] and Table [7] suggest that the IDV properties depend not only on the 5 GHz flux 
density, but also on the optical spectral type. Although the statistical tests indicate a clear differ- 
ence between the distribution of D(2days) in the Type and Type 1 sources for the full sample, 
there are interdependencies that could cause an apparent correlation. 

To study the apparent correlation between spectral type and ISS, we first ask if the different 
redshift distributions of Type and Type 1 sources could cause the difference in the L>(2days). 
From previous st udies, BL Lacs and FSRQs have been found to have different observed redshift 



distributions (e.g. iGiommi et al.ll2012l ) with BL Lacs having a lower median redshift than FSRQs. 
As many (25 out of 59) of our Type sources are without known redshifts, it is probable that 
the redshift distributions of the Type and Type 1 sources are different. Hence it is difficult to 
compare Type and Type 1 sources accurately. 

First we compared the D{2d&ys) distributions of Type sources with known redshift [z^ax ~ 
1.15, N=16) to Type 1 sources with z <1.15 (N=57). The median D(2days) is higher (0.00127 
vs 0.00041) for the radio strong Type sources than Type 1. The statistical tests indicate some 
difference between the samples; pxs = 0.4% that the distributions are similar, and pkw = 2.5% 
that the medians are the same. When the Type 1 redshift cutoff is increased, the difference between 
the two samples increases. No comparisons were made using sources only from the SDSS sample, 
as there are only nine Type sources; however the median D(2days) of Type sources is higher 
than that of Type 1 sources. 

Next we compare all Type radio strong sources (N=31, median L>(2days)=0. 00070) to radio 
strong Type 1 sources with different redshift cutoffs. The KS test indicates p = 2.3% probability 
that the Type 1 with z <1.15, and Type L'(2days) distributions are drawn from the same parent 
population and the KW probability that the medians are the same is pkw = 2.3%. When the Type 
1 redshift cutoff is increased, the KW and KS probabilities that the medians or distributions are 
similar decreases. For example, comparing all Type radio strong sources to radio strong Type 1 
sources with z <1.3 the Pkw = 1.1% and pKS = 0.5%. Considering only radio strong sources from 
the SDSS sub-sample, the median L'(2days) of the Type is also higher than that of the Type 1 
sources, but the difference is not statistically significant for these smaller samples. 

For the radio weak samples, including all Type and Type 1 sources, there is a statistically 
significant difference between the two samples {pks = 0.3% pxw = 1-8%). However, including 
only z < 1.5 Type 1 sources the distributions are similar. Finally, excluding the low redshift [z < 
0.16) Type sources, mainly PEG, and including only z < 1.5 Type 1 sources the statistical tests 
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indicate significant difference between the two samples {pks = 0.3% pkw = 2.1%). 



Using the current data, we conclude that the Type sources have higher Z)(2days) values than 
Type 1 sources. However, the statistical significance is dependent on the sample selection, especially 
the redshift cutoff of Type 1 sources. As many Type sources are without known redshift, it is 
impossible to draw solid conclusions regarding whether there is a real D(2days) difference between 
the Type and Type 1 source s. In general, the red shift distributions of BLL (Type 0) and FSRQ 



(Type 1) are d ifferent (see e.g . iGiommi et al.ll2012l ). with BL Lacs having a lower median redshift 



than FSRQs. iGiommi et al.l (j2012l ) proposed that BLL and FSRQ would have similar intrinsic 
redshift distributions, however, in this case the Z)(2days) distributions of Type and Type 1 do 
not match, suggesting different Z)(2days) properties between Type and Type 1 sources. 

Table [8] shows the number of sources of each spectral type compared with the number of 
M ASIV Survey epoch s in which the source varied according to the visual classification presented 
bv iLovell et al.l (|2008l ). This suggests that Type sources show ISS more often than the Type 1 
sources, which is consistent with the Z)(2days) results above. It is notable that Type sources have 
the highest levels of ISS in both the weak and strong groups, which implies that they typically have 
more compact radio structure than Type 1 sources. The ISS intermittency is likely to be source 
related, rather than purely due to ISM intermittency. These results are relevant to attempts to 
understand the physical differences between the two groups. 



5.2. ISS versus Redshift 



Lovell et al.l (120081 ) found a redshift dependence in the ISS properties. The original 275 radio- 
selected MASIV sources with redshifts revealed that ISS of compact sources at 5 GHz decreases 
with redshift above ~2. From the optical data presented here we can examine whether this decrease 
depends on spectral classification and re-examine the decrease with redshift using the larger sample 
of 320 objects with redshifts. As the D(2days) correlates with the 5 GHz flux density, if the 
5 GHz fiux density were to correlate with redshift, this could cause an artificial -D(2days)-redshift 
correlation. We ran Pearson and Spearman correlation coefficient tests for the full and 'SDSS 
sample' radio strong and weak Type and Type 1 sources. Out of eight correlation coefficients 
only one was statistically significant, the full radio strong sample Pearson test (r ~ -0.18; p ~3%). 
All the other tests indicated no correlation between S^ghz and redshift. 

The 2008 MASIV analysis of Z)(2days) versus redshift included all sources regardless of their 
spectral classification, specifically including Type sources with known low redshift and excluding 
the Type sources with unknown, possibly high, redshift. We ask whether the presence of Type 
sources could be partly responsible for the drop in ISS with increasing redshift. We found (see 
above) the Type sources to be among the strongest scintillators and also to be concentrated at 
redshifts less than 0.7, however almost 50% (27) of the sources do not have a redshift. Fortunately, 
the redshifts of the Type sources are less than 1.15, so the possibly misclassifled Type 1 sources 



will most likely affect the low-redshift end of the Type 1 sources. iLovell et alj (|2008l ) show that the 
average -D(2days) values drop around z ~ 2, hence the misclassified Type 1 sources should have 
negligible bias to the ISS-redshift correlation. 



We now restrict the analysis to the Type 1 objects which constitute the great majority of 
our sources. Figure [8] shows the D(2days) against redshift, with horizontal lines indicating the 
median L)(2days) for different redshifts bins. Both radio strong and weak sources have the lowest 
D(2days) median at the highest redshift bin (z > 2.5). Combining the radio strong and weak 
samples, the KW test indicates 0.5% probability that the 1.5 < z < 2.5 and z > 2.5 samples have 
similar Z)(2days) median and the KS-probability that the two samples are drawn from the same 
parent population is 3.2%. We also note that the two highest redshift bins (D(2days) = 0.00029 and 
0.00051) have about equal numbers of radio strong and radio weak sources, reducing the probability 
that the Z)(2days) decrease with redshift is induced by the 5 GHz flux density. Figure [8] suggests 
that the radio weak sample has a stronger redshift dependence than the radio strong sample. The 
KW probability that the weak sample -D(2days) medians at z > 2.5 and 1.5 < z < 2.5 are the same 
is 2.0% and the KS for the same parent population is 1.4%. From these results we conclude that 
Type 1 sources show a decrease in Z)(2days) with increasing redshift, though the trend is not as 
strong statistically as the decrease versus 5 GHz flux density. 



5.3. ISS and Radio Power 

We searched for correlations between radio power and the -D(2days), and compared the radio 
powers of ISS and non-ISS sources. We calculated the rest-frame radio power at 5 GHz for isotropic 
emission as 

Pr = 47rDlS,{l + zr^'+»^\ (4) 

where Dl is the luminosity distance in the adopted ACDM cosmology, is the mean first epoch 
MASIV VLA flux density and the last term is for the ^-correction. We have assumed ar=0 for all 
our sources. 

As the initial sample selection has a gap in the 5 GHz flux density distribution between about 
0.2-0.5 Jy, the radio power distribution is not continuous at a given redshift. In order to study 
continuous distributions, we selected Type 1 sources with 1.0 < z < 1.8, 59 radio strong and 39 
radio weak sources, and divided the sample into ISS sources (with i:>(2days) > 4 x 10""^) and non- 
ISS sources (L»(2days) < 4x10"^) to compare the radio power distributions. The KS and KW tests 
suggest significant difference in the radio power distributions {pks = 0.9%, pxw = 0.3%), non-ISS 
sources being more powerful. The Pearson correlation coefficient indicates weak correlation (r=- 
0.28 p=0.5%) between radio power and Z)(2days). However, as there are fewer radio weak sources 
than radio strong sources (Table O Figure [6]), this can cause a correlation between ISS and radio 
power. Thus we divided the sample into radio strong and radio weak and compared the radio power 
distributions of the ISS and non-ISS sources. The KS and KW tests show no difference between 
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the samples, suggesting that the correlation seen by the full sample is weak or non-existent. The 
redshift distributions of the samples are identical in each case and the results are not sensitive to 
the selected ISS/non-ISS limit. Our interpretation of the results is that the ISS-P,- correlation is 
weak. 



Finally, we note that using the radio loudness limit of Pr= 10 W/Hz, (similar to lLaurent-Muehleisen et al 



19981 ) five low redshift [z <0.05) sources are "radio quiet 



5.4. ISS and Optical Brightness 

In the optical, we searched for correlations between L'(2days) and optical magnitude/optical 
luminosity. We selected subsamples, as above, based on the radio power and carried out the 
Spearman and Pearson correlation tests. Apart from Type radio weak sources we found no 
evidence of correlation between the optical /^-magnitude and -D(2days). As several Type sources 
have low redshifts, and might have a strong host galaxy contribution to the optical magnitude, we 
excluded the low redshift sources (0.01 < z < 0.15), and the indication of correlation disappeared. 
We calculated the rest-frame optical power as in equation H] but used the R-band flux density and a 
/c-correction — 2.5(1+q:o)^05io(1 + ^;), where we assumed ao=— 0.5 for all our sources. The statistical 
tests (KS, KW, Spearman) indicate no difference in the optical luminosities of the ISS and non-ISS 
sources nor any correlation between the -D(2days) and optical luminosity. 



5.5. Modeling the ISS dependence on Redshift 



The 275 radio-selected MASIV sources with redshifts known prior to the present paper reve aled 
that ISS of compact sources at 5 GHz decreases with redshift above z ~ 2 (ILovell et al.ll2008l ). In 
Figure 3 of lKoay et al.l (|2012l ) we considered the possible influence of a cutoff in radio power or in 
redshift and found that either a single cutoff beyond redshift z=2 explains the data, or different 
cutoffs in Pj. would be required for the weak and strong samples. Here we follow the simpler 
hypothesis of a physical model in which Z)(2days) is governed by S^ghz and redshift r egardless of 



Pr. T he model is motivated by the physics of jet emission as discussed below and also bv lKoav et al 



(|2012l ). This model also explains why the radio weak sample ha s a stronger redsh ift dependence of 
Z)(2days) compared to the radio strong sample (see Figure 9 in lKoav et al.ll2012l ). 



In Figure[9]we have averaged Z?(2days) into four redshift bins for the FSRQs alone, combining 
both weak and strong sources. We note that excluding the Type 2 sources increases the level of ISS 
in the lowest redshift bin and excluding the Type sou rces decreases the level of ISS in the lowest 
redshift bin, when compared with our earlier results (jLovell et al.ll2008l ). We also over plot the 



"'J0248+0434, J1141+5953 and J1604+1744 are classified as galaxies, J1112+3527 as PEG and J0057+3021 as a 
LINER 
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data with curves from a simple theoretical model. The model makes many assumptions and so is 
quite simplistic; however, it has the merit of including the most obvious cosmologically important 
parameters. We assume that the scattering is due to a single layer of isotropic Kolmogorov turbu- 
lence in the interstellar electron density centered at 500 pc from the Earth. We further assume that 
the scintillation for a point source at 5 GHz is at the transition from weak to strong scintillation. 
We model the compact part of the sources by a relativistic jet pointed toward the Earth emitting 
beamed radiation characterized by a brightness temperature in the jet frame of T^o = 10^^ K which 
delivers an observable flux density of 1 Jy. This provides a crude model for synchrotron radiation 
that is self-absorbed causing an upper limit in brightness. The jet is characterized by a Doppler 
factor r which is assumed independent of z. The apparent brightness temperature of the source at 
redshift z is thus given by Tfeor/(l + z). 

Since for a flux density limited sample the apparent angular diameter varies as the inverse 
square root of the apparent brightness temperature, the (1+^;) reduction factor causes the appar- 
ent diameter of the compact source to increase as \/\ + 2, which corresponds to the cosmological 
"angular size versus redshift" phenomenon. VLBI observers have long sought to detect this pre- 
dicted increase in a ngular diameter with redshift, but with only partial success. For example, 
Gurvits et al.l (jl999l ) show angular diameters of compact radio sources that exhibit a barely sig- 
nificant increase beyond a redshift of about 1. We see a highly significant drop in Z)(2days) with 
redshift, but a model for the ISS phenomenon is required before it can be interpreted cosmo- 
log ically. An ISS mod el was applied to the data from the Green Bank monitoring observations 
by iRickett et al.l (120061). However , here we have improved the model by using the formulae de- 
veloped by iGoodman &: NarayanI (|2006l ). which explicitly include contributions from diffractive 
scintillations, th at are important ne ar the transition from weak to strong scintillation and were 
not included by iRickett et al.l (120061 ) . An understanding of D(2days) is obtained by noting its 
asymptotic dependence oc Q^omce fo'^ intraday ISS. This assumes that the ISS is quenched and that 
its time scale is not much longer than 2 days. However, if the time scales are much longer, then 
Z)(2days) drops even more steeply with ^source- The values chosen for T^g, T and the mean flux 
density Sc in the compact source enter through a single quantity Ty^ jSc- Thus, while the model 
exhibits the expected dependence on redshift, we should also expect a wide range in this unknown 
quantity. 

A single model curve (T = 20) can be seen to overlap three of the four observed values within 
±1(7. However, in the highest redshift bin, L'(2days) lies more than 2a below the curve. A steeper 
than predicted drop in ISS with redshift could imply redshift evolution in the Doppler factor of the 
jet s or angular broade ning caused by propagation through the intergalactic medium, as mentioned 
bv lLovell et al.l (|2008l ). However, the discrepancy shown here is not statistically strong enough and 
we await further redshift values from the remaining (weak sample) sources before further analysis. 
In fact, this discrepancy at the highest redshift bin could be linked to a steepening of the mean 
source spectral index with increasing redshift, coupled with the lower Z)(2d ays) of steeper spec trum 
sources, as has been found in the dual-frequency follow-up observations (jKoay et al.ll2012l ). As 
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noted in Section 2, the flux densities used to estimate the spectral indices for the selection of 
MASIV sources were non-simultaneous; and the dual-frequency follow-up observations provided 
more accurate estimates. Indeed any more detailed analysis must also consider the distributions in 
the parameters of the model rather than the single values adopted here. Nevertheless, we stress that 
the present work presents the strongest observational evidence to date for the predicted increase in 
angular size due to cosmological expansion. 

6. SUMMARY 

We present the optical and spectroscopic identiflcations of a sample of 475 flat-spectrum radio 
sources which arc unresolved at 5 GHz (FWHM < 50 mas) . About half of these sources also show 
ISS. Our sample is divided into radio strong {S^gHz > 0.3 Jy) and radio weak ( 0.06< S^gHz <0.3 
Jy) categories. For optical identification the SDSS DR5 and GSC2.3 surveys were used. About 
80% of the radio strong and 60% of the radio weak sources have an optical counterpart with R < 
20 magnitude. Using sub-samples defined from the SDSS data indicates that the nearly complete 
radio strong sample ((Ssghz) = 0.9 Jy, N=78) has median r-magnitude ~ 17.9 (rmax ~ 22.0), and 
the 70%-identifled radio weak sample has r ~ 18.9 (rmax ~ 22.9, (-Ssghz) =0.1 Jy, N=95). 

Spectroscopic identiflcations and redshifts are from the literature and from our own observa- 
tions (347 sources). Spectroscopic identiflcation of the radio strong sample is 91% complete and 
the weak sample is 57% complete. 

The key flndings can be summarized as follows: 

1) The radio weak sources have flatter radio-optical spectral index than the radio strong 
sources. 

2) Almost 80% of the MASIV sources are identifled as broad line objects, ~13% as BL Lacs 
and ~7% as narrow line objects or galaxies. The spectroscopic identiflcations are similar for the 
radio strong and weak samples. Our results suggest that selecting for compact radio structure 
favors FSRQs and BL Lacs over other flat spectrum sources. 

3) Of the scintillating sources, 25% are BL Lacs and the rest are almost exclusively broad 
line objects. The fraction of BL Lacs increases with the degree of persistence of ISS and 40% of 
the sources with most persistent ISS (scintillating at all four epochs) are BL Lac objects. This 
is consistent with BL Lacs being more strongly core dominated than FSRQs, and could possibly 
indicate relativistic jets that are pointed closer to the direction of the Earth. We conflrm that 
D(2days) depends on 5 GHz flux density, with weak sources having higher D(2days) values. A 
similar trend is foTind also when broad line or weak line objects arc stTidicd separately. We found 
indications that the ISS properties correlate with emission line equivalent width. 

4) Radio and optical luminosities of the scintillating and non-scintillating broad line sources 
are similar at a given redshift, suggesting intrinsically similar SEDs. 
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5) We find no correlation between the optical luminosity and ISS. We find a weak correlation 
between radio power and ISS. However at given redshift the ISS and non-ISS sources have similar 
radio power. 

6) For broad line objects, we confirm the sharp decrease in the number of ISS sources and in 
the level of their ISS above a redshift ~ 2. The decrease is compared to a simple model for ISS 
of flat spectrum radio sources with maximum brightness temperatures that are Doppler-boosted 
in jets pointing toward the Earth. There is reasonable agreement but more redshifts are needed 
for the weak radio sample. We found strong observational evidence for the predicted increase in 
angular size due to cosmological expansion. 
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Fig. 1. — Apparent SDSS r-magnitude (red dots) and GSC2.3 F-magnitude (gray circles) versus 
separation between the radio and optical positions. 
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Fig. 2. — The histogram of the magnitude difference for objects with both SDSS and GSC identi- 
fication. The bottom panel shows the converted R-magnitude difference. 
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Fig. 3. — 5 GHz flux density vs optical flux. The blue circles, red stars and open circles are 
for variable, non-variable and intermittently variable, respectively. The magnitude upper limits 
(R=20.5) are indicated with blue, red, black triangles for variable, non-variable and intermittently 
variable, respectively. The diagonal lines indicate constant radio-optical spectral index -0.2, -0.45, 
-0.55 and -1.0 from top to bottom. 
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Fig. 4. — Example spectra of an FSRQ (Type 1), NLRG (Type 2) and BL Lac object. In the top 
panel the vertical bars indicate redshifted Lya, Si II, C IVA1549, C IIIA1909 and Mg IIA2798. In the 
middle panel vertical bars indicate redshifted Ne VA3426, O IIA3730, Ne IIIA3870, [O III]AA4958, 
5007 and Ha. In the bottom panel the vertical bars indicate absorption lines at redshift 1.41 with 
Mg IA2852, Mg IIAA2803,2796 and Fe IIA2344. The e-signs identify features due to the Earth's 
atmosphere 




Fig. 5. — Redshift histograms of MASIV subsamples selected by 5 GHz flux density and optical 
spectral type. The BL Lacs with featureless spectrum are in the first bin of Type panels. The 
vertical line in each subplot indicates the median redshift. 
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Fig. 6. — Optical R-band luminosity versus radio power, with FSRQ, BL Lac and AGN marked 
with red circles, blue squares and green stars respectively. The filled symbols show in the top panel 
the location of optically bright sources, in the middle panel radio weak sources and in the bottom 
panel redshift 1.3 < z < 1.8 objects. Note the two outliers (J0248+0434 and J1141+5953) near 
Mr=-26, whose magnitude estimates from the USNO-Bl are too bright. 
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Fig. 7. — D(2days) histograms of MASIV subsamples selected by 5 GHz flux density and optical 
spectral type. The black vertical line at —3.4 indicates the limit where SF-analysis suggests vari- 
ability. In each subplot the red vertical line indicates the median value and the black line the mean 
value. 
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Fig. 8. — D(2days) against redshift for the Type 1 sources. Red triangles indicate radio strong 
sources and blue circles radio weak sources. The horizontal thin and thick lines indicate the median 
Z)(2days) for the given redshift interval for radio strong and weak samples, respectively. Type 
objects are not shown as most of them do not have redshifts. 
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Fig. 9. — The data for all FSRQs averaged into 4 redshift bins. Over plotted are theoretical curves 
for a simple ISS model, discussed in the text, with Doppler factors as indicated in the legend. The 
horizontal line is at 4 xlO~^, which is the threshold for significant variation. 
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Table 1. Summary of the optical identification using GSC2.3 and SDSS DR5. 



5f;//;_ limit 


N UKxiian EF sur\'{^y 


S>0.3Jy 
S<0.3Jy 


193 18.2 36 GSC2.3 F 
167 19.5 79 GSC2.3F 


7' 


^ <RA< 16^ and S < 64 


S>0.3Jy 
S<0.3Jy 


62 18.25 4 SDSSr 
75 19.82 11 SDSSr 


Note. — 


Columns:(l) the radio flux density 



limit, (2) the number of detected optical counter 
parts (3) the median magnitude of the sample (no 
galactic extinction correction) (4) non-detections 
(5) The survey and pass band 
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COMMENTS ON ARCHIVAL REDSHIFTS 



After careful inspection three objects with archival redshift were rejected. In addition, un- 
reliable redshifts have been proposed for ten BL Lac objects and four objects with adopted, but 
uncertain redshift or spectroscopic identification. We briefly discuss these here. 



J0200+0322: iFricke et al.1 (|l983l ) reported redshift, z=0.765, based on a single (Mg IIA2798) hue 
at 4920A. Their spectrum had the wavelength range from 4200-7100A where for a typical QSO 
spectrum a si ngle broad lin e can be detected at z ~ 0.8 (Mg II) or 1.6 (C III]). 
J0204+1514: Stickel et al.l (|l996l ) reported -2=0. 833 based on two narrow emission lines. Later 



Perlman et al. 



(| 19981 ) found more lines and a better fit for z=0.405 and they classified this object 
as an FSRQ, which we adopted for this work. However this object should be re-observed in order 
to secure the spectroscop i c ideri tification. 

J0406+2511: iHook et all (|l996l ) assumed the single line at 6260A to be [O II]A3727, at z=0.68. 
Their data has a wavelength range from 5200-9100A, however the red and blue ends are very noisy, 
reducing the useful range to 5800-8500A. If the redshift is indeed 0.68, H/3 should be around 8150A, 
and/or [O III ] A5007 around 8410A, whic h are not seen. 



J0449+1121: Ivon Montignv et al.l (|l995l ) reported 2;=1.27, however iHalpern et al.l (|2003l ) found fea- 
tureless spectrum. Thus we adopted a BL Lac identification for this source. 

J0738+1742: NED re ports z = 0.424. Howev er, this is the redshift of intervening Mgll absorption 



feature, i.e. z > 0.43 (jRector &: Stockd (120011) and refere r ices t herein). 

J0745+1011: We ad opted the redshift from iBest et al.l (|2003l ) who measured z=2.624, however, 



Labiano et al.l ( 20071 ) could not confirm this. If the redshift is correct the radio power is 4x lO^^W/Hz 
and an outlier in the M/{-Log(P)-plot (Figured]). This object should be re-observed. 
J0753+5352: This sour ce has featureless optic al spectrum and the redshift has a lower limit of 0.2 
based on imaging data (jStickel &: Kuehrlll993l ). 

J0818+42 22: This is a rad i o sour ce with featureless optical spectrum and has no reliable redshift 
available (jRector Stockd (j200ll ) and references therein). 

J0929+5013: The SDSS redshift is uncertain as there are n o obvious emissio n lines . 
J1150+2417: This object has a featureless optical spectrum iRector &: Stockd (|200ll . and references 
therein). 



J1309+1154 This object has a featureless optical spectrum by iLvnds Willa (jl972l ) and SDSS. 
Note that this object has SDSS DR5 redshift of 0.3852 a nd SDSS DR 8 redshift of 2.60 1. 
J1446+1721: The redshift of this FS RQ is uncertain as ISowards-Emmerd et al.l (|2005l ) reported, 
2;=0.102, however iHealey et al.l ([200a) reported 2;=1.026. For this work we adopt FSRQ classifica- 
tion with unknown redshift. 



J1502+3350: IWhite et al.l (I2OOOI ) report featureless spectrum for this object. Note that the SDSS 
DR6 suggest 2=2.222, however the spectrum appears fe atureless. 



J1602+3326: Th e redshift is based on opt ical magnitude (jSnellen et al.ll200d ) using the GPS galaxy 
Hubble diagram (jStanghellini et al.ll2005l ). 

J1648+2141 The SDSS spectrum is featureless, however the SDSS DR6 suggest z=1.085. 



Table 2. MASIV Data 



Source 


RA 


Dec 


SsGHz 


Var 


D(2days) 


E{B-V) 


R 


Ref 


z 


Type 


ID 


Ref Comment 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 




J0005+3820 


00:05:57.1755 


38:20:15.168 


0.61 


1 


0.0897 


0.08 


17.00 


G 


0.229 


Tl 


fsrq 


1 


JOOlO+1724 


00:10:33.9920 


17:24:18.791 


0.84 





0.0378 


0.03 


16.96 


G 


1.601 


Tl 


fsrq 


2 


J0017+5312 


00:17:51.7596 


53:12:19.125 


0.90 


2 


0.4006 


0.18 


17.10 


G 


2.574 


Tl 


fsrq 


3,4 


J0017+8135 


00:17:08.4750 


81:35:08.137 


0.61 





0.1857 


0.40 


15.49 


G 


3.387 


Tl 


fsrq 


5 


J0019+2021 


00:19:37.8545 


20:21:45.570 


1.15 


1 


0.0604 


0.06 


18.69 


G 




TO 


bll 


6,7 


J0019+7327 


00:19:45.7871 


73:27:30.019 


0.44 


1 


0.3216 


0.32 


17.48 


G 


1.781 


Tl 


fsrq 


8 



Note. — Columns are as follows (1) lAU name (J2000.0), (2)R.A.(J2000.0), (3)Decl.(J2000.0) (4) first epoch MASIV 5GHz flux density, (5) 
The number of epochs of variability based on visual classification, where 'C indicates secondary calibrators, (6) Structure function _D(2days), (7) 
Galactic extinction, (8) R-magnitude, (9) Reference; G=GSC2.3, S=SDSS DR5, U=USNO-B R2 (10) Redshift, (11) Spectral type (12) Spectroscopic 
identification, (13) Reference for redshift and spectroscopic identifica tion, (14) Comment 

Re ferences: (l) lstic kcl & Kuchr (1994) (2) Wills fc W ills (19761 (3) ISa,rgent et al.l {HH), (4) iKuhr et"ai] il98i ). (5) ISowards-Emmerd et al.l (|2005l '), 
('6) IChu et aU l ll986l) . (7) Palomar 2007 this paper, (81 [Lawrence et al.l l ll98^ 

Table 1 is published in its entirety in the electronic edition of n, A portion is shown here for guidance regarding its form and content. 
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Table 3. Objects with single emission line from our data (that have an ambiguous 
identification). The range indicates the useful wavelength range of the spectrogram. 



Source 


Line 


Range 






A 


A 


per pixel 


J0213+3652 


4179 


3500-5000 


2 


J0342+3859 


5453 


4150-7500 


15 


J1316+6927 


5675 


4500-7500 


3 


J1711+6853 


5660 


3900-8500 


8 


J1953+3537 


5125 


5200-7500 


10 


J2012+5308 


3967 


3950-7500 


6 



Table 4. Summary of the spectroscopic identification. The third and fourth column indicates the 
number of sources in the radio strong and radio weak samples, respectively. 



Type 


Spect id 


Strong 


Weak 


TypeO 


BLL 


37 


11 




BLL? 




5 




PEG 




6 


Type 1 


FSRQ 


155 


107 




FSRQ? 


1 






Sy 1 




2 




Sy 1/BLRG 


1 






NLSyl 


1 






NLQSO 


2 






Syl.5/LINER 


1 




Type 2 


Sy2 


3 


2 




Sy2? 




1 




NLRG 


5 


1 




LINER 


1 


1 




HIT? 




1 




galaxy 




3 


Single line objects 




1 


5 
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J1649+7442: This redshift is from lAppenzeller et al.l (|l998l ) from the spectrum with a wavelength 
range from 4000 to 7000/9000 A. The redshift is based on a single line at 6658 A which is assumed 
to be Mg IIA2798 at ^=1.378. At this redshift, one would expect to see C III]A1909 at 4532A, 
however if the line was C III] or H/3 one would expect to see C IV at ~5392A or [O III] near 6870A. 
The redshift i s adopted, but needs to be con firmed. 



J1719+1745: ISowards-Emmerd et alJ (|2005l ). reported z=0.173. However IShaw et al.l (|2009l ) find 
featureless spectrum and they estimate redshift to be >0.58. This source is classified as a BL Lac 
object without known redshift. 



B. TRANSFORMATION FROM GSC2.3 TO LANDOLT SYSTEM 



In order to transform the GSF2.3 magnitudes to the iLandold (|2009l ) system, we obtained 
J FN magnitudes from the Landolt standard stars. We included only stars with four or more 
measurements by Landolt, then plotted J — B, F — R and N — I and removed outliers. Least 
squares polynomial fit resulted following transformations: 



5g5C7 = J + 0.18(J-F)-0.08, 



(Bl) 



Rqsc = F[-0.085(F - A^) + 0.04], 



(B2) 



Icsc = A^[+0.066(F -N)- 0.02]. 



(B3) 



The terms in square brackets for Rgsc and lose a-re from the best least-squares fit, however 
the fit is improved only marginally, hence we have omitted the color term. The root mean squared 
error for the BRI-conversions are 0.10, 0.11 and 0.11, respectively. 

The Base transformation is valid for the (J — F)-colors from -0.5 to 2.5 and Rgsc-, Igsc for 
the {F — A^)-color range from -0.4 to 1.4. 



C. NOTES ON INDIVIDUAL OBJECTS 



In this appendix, we discuss some of our spectroscopic follow-up targets. For the narrow line 
objects we report the EW of the strongest line. 



J0019+2021 BL Lac classification is from IChu et al.l (|1986l ) who found no strong emission lines. 
Our DBSP data (S/N<35) shows absorption system at z=lAl with Mg IA2852, Mg IIAA2803,2796, 
(Fe IIA2344) and a possible absorption system at z=2.274 (Figured]). 
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J0150+2646 The DBSP (07-2007) and NOT (8-2004) spectrum are featureless with S/N<14 and 
S/N<8, resp ectively. Thus, this s ource is tentatively identified as a BL Lac object. 
J0509+0541 iHalpern et al.l (j2003l ) reported featureless spectrum from several observations. Our 
high S/N (>130) data is featureless showing only galactic absorption lines and unidentified absorp- 
tion at 42761. 

J0800+4854 Our data show strong absorption at the C IVA1549- and Lya-lines and also the 
He IlA1640-line is easily detected. 

J1118+2922 This low redshift object has two easily detected narrow emission hues at A3998A and 
A7038A, identified as [O II]A3727 and Ha at redshift of 0.072. The line ratios suggest HH region 
or Liner classification. 

J1352+3603 We detect narrow [O III]AA4958, 5007 and [O II]A3727 lines at z=0.3057 with, Mb=- 
21.5, (Mij=-23.0), C=42%, EW([0 II])=20A and noisy EW(Ha)~ 37A. The source is classified as 
P EG based on the C a-break and line EW. On the other hand, based on the classification scheme 
of lLandt et al.l (j2004l ). this is a strong- lined source. 

J1442+0625 Our data show a single hue at 4753A, with continuum from 4800-7500A (S/N<10). 
Our tentative identification is Mg IlA2798-line at 2;=0.698. If this line would be C III] at z=1.49 
we would expect to see Mg IIA2798 at about 7000A which is not detected. 

J1444+0257 This object has a featureless spectrum from 4500 to 7500A with S/N<5, possibly a 
BL Lac object 

J1444+1632 We detect a sing le line at 5688A with noisy continuum (S/N <3) and tentative line 

at 7247A, these are identified as Lya and C IVA1549, which would indicate z=3 .7. 

J1505+0326 This object has been classified as FSRQ in CGRaBS (jHealev et al.lboosl 'l. Our NOT 
data confirms the CGRaBS and SDSS redshift. However reanalyzing the SDSS data and based on 
the NOT data the rest frame line widths are less than 2000km/s indicating NLSyl classification. 
J1651+0129 Our DBSP data shows a possible line at 7128A with the continuum S/N<7. As we 
are not able to identify the line or spectrum, this target is not included to the analysis. 
J1718+4448 The DBSP red spectrum show no obvious lines from 5300 to 8300A with S/N < 5, 
possibly a BL Lac object. 

J1728+1931 We detect narrow ([O II]A3727), [O III]AA4958, 5007 and a noisy [O I]A6300 lines at 
z=0.1759, with [O IIl]A5007/H/3>3, EW([0 II])=13A, Mb=-21.1, {Mr=-24.2) and C=42%. The 
Ha-line is noisy, howev er it appears to be narrow. This is border line object between strong and 
weak lined sources using lLandt et al.l (|2004l ) classification. We tentatively classify this as a Sy2-type 
object. 

J1734+3857 lStickel et al.l (|l989l ) classification for this source is an FSRQ at z=0.97 based on a sin- 
gle broad emission line which was assumed to be Mg IIA2798 (EW(Mg II)=15.3A, observed) Our 
data shows also only a sing le line at 5535A assumed to be Mg II at z=0.97 with the EW=4.2A, 

hence BL Lac classification. 

J1739+4737 This is a BL Lac-object with featureless spectrum ([Scarpa et al.ll2000l ). Our DBSP 
data show no lines with S/N<20 

J1742-I-5945 This object has a featureless spectrum iLaurent-Muehleisen et al.l (|l998l ). and the red- 
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shift {z=0.4:) is based on a detection of the host galaxy (jNilsson et al.l 120031 ) . Our NOT spectrum 
(07-2006) shows no hues from 4000 to 7500 A with S/ N<25. 

J1747+4658 This source has no emission Hne redshift. IVermeulen et alj (|l996l ) reported an inter- 
vening Mg II absorption line system at z= 1.484, which, by some authors, has been adopted as the 
source redshift. Our NOT spectrum (Grism 6) from 4000 to 5500A with S/N<6 is featureless. 
J1757+0531 Narrow emission lines, [O II]A3727 and [O III]AA4958, 5007 at z=0. 34482 are clearly 
visible and so is a noisy strong Ha feature at the red end of the spectrogram. The line and continuum 
properties suggest Type 2 identification: [O IIl]A5007/H/3>3, Mr=-23.2, EW([0 III]A5007)=37A and 
C=28%. However the continuum is rather noisy, (S/N<8) and weak broad emission lines might 
have been undete cted (Figure [H) . Tentati ve spectroscopic identification is NLRG, however, based 
on Fig. 8 and 9 of lCaccianiga et al.l (|2008l ) this object could be a Type 1 object with the AGN core 
diluted by the host galaxy. 

J1832+1357 Our data confirms the CGRaBS (jHealev et al.l l2008l ) redshift. We detect also some 
absorption f eatures in C IVA 1549-line 

J1848+3219 ISowards-Emmerd et al.l poosl ) reported this source as an FSRQ at 2;=0.798. Our 
data shows a strong line at 5035 A, a noisy narrow hne at 6703A (identified as Mg IIA2798 and 
[O II]A3727, respectively) and no strong hues in 3600-5000 and 5100-7500A with S/N<17, sup- 
porting the literature identification and redshift. 

J1905+1943 A strong broad line is detected at 5106A and a tentative broad line at 6267A, sug- 
gesting z =2.3. 

J1950+0807 From the DBSP data we detect narrow [O II]A3727, [O III]AA4958, 5007 and broad 
Ha-lines at ^=0.2964, with Mr=-22.8, Mb=-20.7 and EW([0 III]A5007)=-48A, suggesting Syl.9 
classification. 

J2023+5427 The redshift 1.48 is measured from Mg IIA2798 and C IIl]A1909-lines. In addition 
we detect Mgl(2852), MgII(2803,2796), FeII(2387,2600) and FeII(2383,2344) absorption lines at 
z=1.415 

J2130+0339 From the DBSP data we detect narrow hues [O II]A3727 and [O III]AA4958, 5007 at 
z=0.653, with EW([0 III]A5007)=140A, Mr >-22.5 and C=31%, suggesting Type 2 identification. 
Note, this source is not detected in the GSC2.3. 

J2155+0916 This object is tentatively identified as a BL Lac object. Our NOT spectrum is fea- 
tureless from 4500 to 7500A with S/N<20, however if the redshift is ~0.6, both Mg IIA2798 and 
H/3-lines will be outside the spectrogram. 

J2201+5048 We confirm the lSowards-Emmerd et al.l (|2005l ) redshift, 1.899 and in addition we detect 
Mg Il-absorption system at z=1.3822 with Mg IA2852 MgIlAA2803,2796, Fe IIA2344,2374,2383, 
however this source was rejected from the iLovell et al.l (120081 ) analysis due to evidence of source 
structure, h ence not included t o the analysis. 

,12203+1725 ISmith et aP (Il977l ) classified this source as a FSRQ based on strong C III]A1909- and 
Mg IlA2798-nnes at z=1.076. Our data shows a single line at 581lA with EW(Mg II)=3.3A and 
S/N<60 around the line, hence BL Lac object classification. 

J2212+2759 This is a BL Lac object, with no detected emission lines, based on our DBSP data 
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(07-2007) with S/N<45 and NOT data (07-2005 ) with S/N<18. 
J2230+6946 We confirm the lHealev etaP ^20od ) (CGRaBS) redshift. 

J2237+4216 An emission hne is detected at 5655A and this is the only emission hne between 4500 
and 7500A, suggesting Mg II ident ification. 

J2241+4120 iHenstock et alj (jl997|) and ISowards-Emmerd et al.l ((20051) classified this as a BL Lac 
object. Our DBSP data (07-2007) with S/N<60 and NOT data (07-2005) from 4000 to 7500A with 
S/N<30 confirm this. 

J2242+2955 A strong broad line is detected at 4102A and a tentative broad line at 5013A, sug- 
gesting z =1.7. 

J2258+0516 We detected a strong absorption feature at the centre of the C IVA1549-line and Lya 
line is almost fully absorbed. 

J2303+1431 This is a BL Lac object with featureless spectrum. Our DBSP spectrum (07-2007) 
has S/N<25 and NOT spectrum (07-2006) S/N<8. 

J2315+8631 A broad hue at 5273A (Mg IIA2798) is visible in our DBSP (08-2007) and NOT 
(07-2006) spectrum. This is the only strong emission line between 3700 and 8500A (S/N<30) , 
hence a tent ative redshif t of z= 0.88. 

J2325+3957 1 Shaw et al.l (120091 ) reported a featureless spectrum with intervening Mg IIA2798 and 
Fell absorption features at z=1.04. Our DBSP data (07-2007) shows no emission lines. However, 
we detect absorption features which are identified as Mg IA2852, Mg IIAA2803,2796, Fe IIA2600 
(and Fe IIA2344) at z=0.4695 and Mg IIAA2803,2796, Fe IIA2600 at z=0.4156 with the continuum 
S/N up to ~25. We also detected an absorption feature at 5757A possibly Mg IIA2798 at z=1.04, 
but not the expected Fe II feature at 5353A. The 4119A feature, Mg IIA2803 at z=0.4695, is also 
visible from the NOT spectrum (07-2006). 
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Table 5. Classification break down of the radio strong and weak samples including the secondary 
cahbrators. The 'SDSS sample' source distributions are in column four and five. 



type 


^full % ^SDSS 


% 




Radio strong sources 


Type 
Type 1 
Type 2 


37 18±3% 19 
162 78±6% 57 
9 4±1% 2 


24±6 % 
72±10% 
3±2 % 


Radio weak sample 


Type 
Type 1 
Type 2 


22 16±3% 8 
108 78±7% 46 
9 6±2% 5 


14±5% 
78±11% 
8±4% 



Note. — The errors are Poissonian. 
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Table 6. Summary of some flat radio spectrum AGN surveys. Note, the PEGs in CBS and 
MASIV are included in the "other AGN" as is the case with CGRaBS and DXRBS. 



Survey 


5 GHz hmit 


R-mag 


FSRQ 


BL Lac 


other AGN 


Comment 




mJy 




% 


% 


% 




CBS 


30 


17.5 


45 


21 


34 




CBS rb 


65: 


17.5 


51 


20 


29 


MASIV 5 GHz hmit 


CGRaBS 


65 


23.7 


88 


10 


2 


Selection by Figure of Merit 


DXRBS 


50 


24.0 


80 


12 


8 


some steep spectrum sources 


MASIV 


60 


23.1 


78 


15 


7 


Selection by 8.4 GHz 


MASIV ob 


60 


17.5: 


64 


24 


12 


Optically bright CBS limit 


MASIV ISS 


60 


21.3 


64 


34 


2 


Persistent IDV sources, IDV3, 4 


MASIV SDSS 


60 


21.7 


75 


18 


7 


SDSS sample 



Note. — Columns: (1) the name of the sample, (2) 5 GHz radio flux density limit for faintest objects, (3) 
optical magnitude limit or the faintest detected objects, (4) the fraction of broad line objects, (5) the fraction 
of BL Lac objects (6) the fraction of other spectral type objects, mostly narrow line objects. 
The colon indicates MASIV radio flux density limit for the "CBS rb" , and CBS opt ical limit for "MASIV ob" . 
Note, the DXRBS includes steep radio spectrum sources as well (ILandt et al.ll200ll ). 



Table 7. Median Z)(2days) of the subsamples selected by optical spectral type and 5 GHz flux 
density. The First column indicates optical spectral type, followed by the median L>(2days) and 

the number of objects. 



Type 


D(2days)xl000 


N 


L»(2days)xl000 


N 




S5GHZ > O.SJy 




S5GHZ < 0.3Jy 







0.701 


31 


2.35 


22 


1 


0.335 


143 


0.868 


108 


2 


0.193 


6 


0.44 


9 




Type 1 z 


< 1.4 sources 





1 0.346 79 1.27 54 
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Table 8. Number of sources and the fractions of "variable" MASIV epochs for the three optical 
spectral types. The first column indicates the number variable in epochs and the next columns 
show the fraction of sources in each group with the number of sources in parentheses. 



Sample 


TypeO 


Type 1 


Type 2 


IDV 


strong weak 


strong weak 


strong weak 





13% (4) 18% (4) 


36% (51) 36% (39) 


5 7 


1 


19% (6) 18% (4) 


24% (34) 15%(16) 


1 1 


2 


13% (4) (0) 


24% (34) 17%(18) 





3 


29% (9) 18% (4) 


8% (12) 18%(19) 


1 


4 


26% (8) 46% (10) 


8% (12) 15%(16) 


1 



Table 9. Redshift, core radio power and optical luminosity, excluding secondary calibra tors. The 
IDVO, IDVl, etc. indicates the number of epochs of variability (jLovell et al.ll2008l ). 



Strong sources Weak sources 









1\T 


<^ z ^ 






l\ 


Z ^ 




























TypeO 






16 


0.4±0.3 


26.6±0.6 


-25.2±1.1 


10 


0.2±0.1 


25.0±0.8 


-23.6±1.3 


Type 1 






142 


1.3±0.7 


27.6±0.6 


-26.4±1.6 


108 


1.4±0.8 


26.8±0.5 


-25.9±1.6 


Type 2 






6 


0.4±0.3 


26.7±1.3 




9 


0.1±0.2 


24.1±1.1 
















Type 1 










L>(2days)> 4 x 


10" 


-4 


59 


1.2±0.6 


27.5±0.6 


-26.2±1.4 


76 


1.3±0.7 


26.7±0.5 


-25.7±1.5 


L'(2days)< 4 x 


10- 


-4 


83 


1.4±0.7 


27.7±0.5 


-26.7±1.6 


32 


1.8±1.0 


27.0±0.5 


-26.7±1.4 














1.0 < z < 1.8 










L'(2days)> 4 x 


10" 


-4 


24 


1.4±0.2 


27.5±0.2 


-26.4±1.0 


29 


1.3±0.2 


26.7±0.2 


-26. Oil. 2 


Z)(2days)< 4 x 


10- 


-4 


35 


1.4±0.2 


27.6±0.2 


-26.4±1.3 


10 


1.3±0.3 


26.8±0.2 


-25.7±1.0 














Type 1 










ISS (IDV234) 






58 


1.2±0.7 


27.5±0.6 


-26.3±1.4 


53 


1.2±0.6 


26.7±0.5 


-25.4±1.4 


IDVO 






50 


1.4±0.8 


27.7±0.6 


-26.7±1.8 


39 


1.9±1.0 


27.1±0.5 


-26.5±1.6 


IDVl 






34 


1.3±0.6 


27.5±0.5 


-26.3±1.3 


16 


1.4±0.7 


26.8±0.4 


-26. Oil. 3 


IDV2 






34 


1.5±0.7 


27.6±0.6 


-26.4±1.5 


18 


1.3±0.7 


26.7±0.5 


-26. Oil. 3 


IDV3 






12 


1.0±0.5 


27.4±0.6 


-25.9±1.2 


19 


1.1±0.7 


26.6±0.6 


-25.0±1.4 


IDV4 






12 


1.4±0.7 


27.6±0.4 


-26.5±1.3 


16 


1.2±0.5 


26.7±0.4 


-25.3±1.4 



Note. — Columns: (1) the sub-sample, (2) Number of objects, (3) Median redshift, (4) Median of 
the Log radio power with standard deviation, (5) Median absolute magnitude with standard deviation, 
(6)-(9) same as Cols. (2)-(5). 



